Submitted to ApJ 

Preprint typeset using I^T^^ style cmulatcapj v. 03/07/07 



O 

o 

< 



O 

U 

6 



> 

00 

cn 
p 

00 

o 

o\ 
o 



X 



AN IMPROVED [O III] LINE WIDTH TO STELLAR VELOCITY DISPERSION CALIBRATION: CURVATURE, 
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ABSTRACT 

An improved transformation of the full width at half maxima (FWHM) of the [O HI] A5007 line in 
AGNs to the stellar velocity dispersion, cr* , of the host galaxy is given. This significantly reduces the 
systematic errors in using FWHM([OIII]) as a proxy for cr*. AGN black hole masses, M,, estimated 
using the Dibai single epoch spectrum method, are combined with the new estimates of cr* to give a 
revised AGN M, -cr* relationship extending up to high masses. This shows that for the most massive 
black holes, M, is systematically higher than predicted by extrapolation of M, oc a* to high masses. 
This supports recent suggestions that stellar dynamical masses of the most massive black holes have 
been systematically underestimated. The steepening of the M, -cr* relationship is consistent with the 
absence of very high cr* galaxies in the local universe and with the curvature of the Faber- Jackson 
relationship. There appears to be significantly less intrinsic scatter in the M, -cr* relationship for 
galaxies with M, > lO^M©. It is speculated that this is connected with the core elliptical versus 
extra- light elliptical dichotomy. The low scatter in the high end of the M, - cr* relationship implies 
that the transformation proposed here and the Dibai method are good indicators of cr* and Af, 
respectively. There is no evidence for evolution of the M, -tr* relationship over time. 
Subject headings: galaxies: active — galaxies:bulges — galaxies: evolution - galaxies: fundamental 
parameters — galaxies: nuclei — quasars:emission lines 



1. INTRODUCTION 

There is considerable interest in the correlation of the 
masses, M,, of supermassive black holes (SMBHs) with 
the stell ar velocity dispersion, , of the bulge of th e host 
galaxv (iGe bhardt et al."200d Iferrarese &: MerrittI 120001 : 
iKormendv & Gebhardt 2001[) because of its important 
implications for the formation and evolution of galax- 
ies and their SMBHs. It is widely recognized that for 
understanding the origin of the M, - cr* relationship it 
is particularly important to know the masses of black 
holes in the most massive and least massive galaxies, 
and to study the evolution of the M, - relationship 
over cosmic time. Having the correct form of the rela- 
tionship is also important for constructing the black hole 
mass function in the local universe and reconciling this 
with the luminosity function of AGNs, energy-generation 
efficiencies, relative accretion rates (Eddington ratios), 
AGN duty cycles, and the evolution of all these quan- 
tities over the history of the universe. The predictions 
from the bulge luminosity, Lbuige, and from the M, -cr, 
relationship conflict for the highest-luminosity galaxies 
because the Af, - Lbuige relationship predicts that the 
brightest galaxies have Af, ~ 10^^ Mq, while the Af, - 
( Tj, relationship predicts Af, < 3 x lO^Af fri at all times 
(jBernardi et al.ll2007bl : iLauer et al.ll2007t ). The number 
of high-mass SMBHs predicted by the Af^, - a* relation- 
ship is also in conflict with the number of relic high-mass 
SMBHs predict ed by the volume density of the most lu- 
minous AGNs (jLauer et al.ll2007t ). 

By far the easiest way to estimate black hole masses 
is to use type-1 (near face-on) AGNs. To do this 
one needs the velocity dispersion of the broad-line 
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region (BLR) gas, which can be readily obtained 
from the H/3 line width, and an effective radius. It 
is also important to establish that the motions of 
the gas being used are dominated by gravity. The 
most direct way of determining an effective radius is 
via "reverberation ma pping" using light echoes fol- 
lowin g AGN variability ( Lvutyi fc CherepashchukI 119721 : 



ityi 

Cher epashchuk fc Lvutvilll973f ). This is most readily ac- 
complis hed by cross corr elating line and continuum vari- 
ability (jGaskell fc Sparke 19861. Velocity-resolved hght 
echoes have sho wn that BLR gas motions are gravitation- 
ally dominated ('Gaskell"198 8HKoratkar fc Gaskel]|ll989l : 
see iGaskell fc Goosmann 20081 and iGaskelll l2009l for re- 
views of the evidence for this) and thus permit estimation 
of Af, from the BLR. Since then, Af, has been estimated 
from reverberation mapping fo r several dozen AGNs (see 
IVestergaard fc PetersonI [20061 for a recent compilation). 
An external check on the accuracy of reverberation map- 
ping masses is provided by the tightness of the correla- 
tions of the Af, estimate s with g„ and Lbuiae - The bulge 
luminosities obtained bv lBentz et al.l (|2009a[ ) imply that 
the error in log Af, determined by reverberation mapping 
is less than ±0.33 dex. 

Reverberatio n map p ing is unfortunately very labor 
intensive, but iDibail (|l977f ) argued that BLR radii, 
and hence masses, could be inferred indirectly from 
photoionization considerations from single-epoch spec- 
tra. Mass estimates from a large nu mber of spec- 
tra o f NGC 5548 in high and low states (jPennev et al.l 
|2009| ) show that random and systematic measuring er- 
rors in logAf, deter mined by the Dibai method are 
small (f« 0.11 dex). iBochkarev fc Gaskelll (l2009l] have 
shown that Dibai's original mass estimates ()Dibailfl977L 
[19801 [1981 are in good agreement with reverberation- 
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mapping mass estimates. The same is true for 
more recent mass determinati ons by the Dibai method 
(jVestergaard fc Peterson|[2006[ ). As with reverberation- 
mapping-determined masses, an independent external 
check is provided by corre lations with the properti es of 
the host-galaxy bulges. .Gaskell &: Kormendvl ()2009l ) find 
a scatter of ±0.23 in logM,. These studies suggest that, 
with care, the Dibai method, which is about two orders of 
magnitude less labor intensive than reverberation map- 
ping, gives M, accurate to < ±0.25 dex. 

A problem in investigating the Af, - a, relationship 
for AGNs is that strong optical continuum continuum 
emission from the AGN accretion disk can make mea- 
suring (T, difficult. This is exacerbated as one goes to 
higher redshifts as the most readily observable stellar 
absorption lines move into the IR. N elson (,20001 pro- 
posed using FWHM[o///]/2.35 as a proxy for be- 
cause the [O III] lines are strong and readily observable. 
For a Gaussian velocity dispersion, a =FWHM/2.35. 
By using FWHM[o ///j/2.35 Nelson showed that active 
galaxies also obey the Af, ~ct» relationship. The use of 
FWHM[o///] as a proxy for cr, has subsequently been 
widely adopted because the [QUI] lines are readily ob- 
servable and they allow studies of the - cr, relation- 
ship to be extended to high er redshifts and higher lumi- 
nosities (jShields et al.ll2003D . 

An obvious concern with using FWHMjq///] as a 
proxy for ct* is that gas, unlike stars, is highly suscepti- 
ble to non-gravitational forces and dissipation, and it has 
long been known that the motions of the gas in galax- 
ies are complicated. In this paper I show that using 
FWHM[o ///] as a proxy for ct* underestimates ct» for 
narrow [O III] lines and overestimates it for broad [O III] 
lines. I give a simple prescription for removing this bias, 
and discuss the implications for the M, -cr» relationship. 

2. USING THE FWHM OF [O III] AS A PROXY FOR 
STELLAR VELOCITY DISPERSION 



The justification 
FWHM[o///]/2.35 as 



of 



iNelsonI (l200nh in using 
a proxy for a* was that in 
the compilation of iNelson fc Whittl j ()1995D and 
iNelson fc Whittle! the average FWHM[n ///]/2.35 
was in good agreement with their averag e directly mea- 
sured a^, . A numbe r of other studi e s (e.g . , iBoroso nl[200l 
'Onken et all 12004 iGreene fc Hoi |2005*; 'Bo nning et all 



2005fc have supported the use of FWHM[o///]/2.35 as 
proxy for a*, although one with significant scatter. 

been known 
that the 



a proxy tor a* 

This scatter is not surprising as it ha s long 
(|Burbidge. Burbidge. fc Prendergas"3 |195£ 
narrow emission lines are blueshifted and this has been 
shown in a number of cases to be the result of outflow. 
The widths of the narrow-line region (NLR) lines also 
depend on th e ionization potential and critical densities 
of the lines (jde Robertis fc Osterbrockl Il986( l. so dif- 
ferent lines will give systematically different FWHMs. 
With these complications we therefore expect a merely 
empirical correlation between FWIIM[(9 777]/2.35 and cr, 
at best. 

Fig. 1 shows FWHM[0 7j7]/2.3 5 and measured cr* val- 
uesfrom lNelson fc Whittii (|1995[ ) and lNelson fc Whittle! 
()1996l ). A dditiona l or u pdated cr* measurements were 
taken from [Nel son (2000) (Mrk 79, Mrk 110, Mrk 279), 



2.6 

2.4 

6 2.2 
o 

-' 2.0 
1.8 
1.6 




• Type-1 AGNs 
OType-2 AGNs 
+ Radio-loud type 1 
X Radio-loud type 2 



Nelson et al. 



Onken et al 



2004 ) (NGC 
20041 ) (NGC 



4151 and NGC 
7469, Mrk 817, 



5548), and 
Ark 120). 



1.6 1.8 2.0 2.2 2.4 2.6 2.8 
Log FWHM III] ^^-35 

Fig. 1. — Measured stellar velocity dispersion, (t», as a function 
of FWHM[o777]/2.35 for AGNs. The dotted line is the equation 
(T* = FWHM[Q 77J] /2.35 and the solid line is an OLS-bisector fit 
lllsobe et aLltloM) for the type-1 AGNs only 

Updated FWHM[o 777] measurements were taken from 
INelson! (|2000l ) for NGC 7469, Mrk 79, Mrk 110, Mrk 279, 
Mrk 817, and Ark 120. Strong rad io sources (objects in 
Table 7 of INelson fc Whittle! Il995l ) are indicated by the 
addition of a "+" to type-1 objects and a "x" to type-2 
objects. 

It can be seen from Fig. 1 that while FWHMfo 777]/2.35 
for all AGNs agrees on average with for the sample as 
a who le, there is substantial scatter. INelson fc Whittlj 
(|1996! ) find an rms scatter of ±0.20 dex (see their Fig. 
7a). The rms scatters for both the type-1 and type-2 
AGNs in Fig. 1 are similar (±0.20 dex ±0.19 dex re- 
spectively). The rms of the quoted measuring errors of 
cr* is ±0.07 dex. Comparison of FWHMs of [QUI] with 
th e FWHMs of simi l ar ion ization [S III] lines (Table 3 
of INelson fc Whittle! 11995! ) suggests that the rms error 
in the measurement of each line is ±0.06 dex. The ex- 
pected scatter in Fig. 1 from measuring errors alone is 
thus expected to be « ±0.10 dex. 

In addition to the large scatter in Fig. 1, it can 
be see that there is a systematic deviation from cr* 
= FWHM[o///]/2.35 for wide [QUI] lines as has al- 
ready been noted by INelson fc Whittle! (|1996l ). It is, 
of course, important for studying the Af, -ct* relation- 
ship for AGNs that there be no system atic errors in 
estimating cr*. INelson fc Whittle! (l996D have already 
pointed out that the systematic deviation is strongest 
for radio-loud AGNs. Fig. 1 shows that it is actually 
the type-2 radio-loud AGNs which show the largest sys- 
tematic effect. This argues that the NLR outflows are 
preferentially near the equatorial plane. Typc-1 AGNs 
show a similar systematic bias when [O HI] is used as 
a proxy for g.^. T his is shown by the OLS-bisector fit 
(|Isobe et al.!ll990! ) to the type-1 AGNs alone in Fig. 1. 
The effect is quite substantial. For type-1 AGNs with 
FWHM[0 777]/2.35 < 100 km s"^ the median under pre- 
diction of (T* is 0.22 dex. Since Af, oc cr^, this translates 
into a 0.9 dex systematic error logAf,. 

The NLR line widths are consistent with other observa- 
tions showing that emission-line gas velocities in bulges 
tend to be sub-virial. For example, in a l arge study 
of elliptical galaxies by Phillip s et al.! (|1986! ) the mean 
FWHM/2.35/cr* from their [N'H] FWHMs is 0.73 ±0.03. 
Spatially resolved Hubble Space Telescope spectroscopy of 
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Fig. 2. — Measured stellar velocity dispersion, (T -jf J clS Qu function 
of the corrected [OIII] line velocity dispersion, c^yo lll]corr^ from 
Eq. 1. The diagonal line is the equation cr* = ct^q j j ^^^^.^ 



NLRs shows that the innermost g as is commgnly m oving 
more slowly than gas further out ()Rice 

Sub-virial NLR motion is obviously being seen because 
gas is settling into a flattened distribution. Since type-1 
AGNs are viewed close to face-on, this will reduce the 
observed FWHM for gas in a flattened distribution. The 
smaller systematic offset for type-2 AGNs would be be- 
cause they are viewed more obliquely. It can be seen in 
Fig. 1 that both the overestimation of ci* and the differ- 
ence between type-1 and type-2 AGNs decreases as the 
FWHM gets wider. The reason for the decreasing differ- 
ence between type-1 and type-2 AGNs with increasing 
FWHM could be that for type-1 AGNs with the largest 
FWHMs the torus opening angle is larger than average, 
and so they can be seen at higher inclinations. The extra 
width of [O III] could also be associated with an outflow. 
This must certainly be going on for the AGNs with the 
highest [OIII] FWHMs. 

3. AN IMPROVED RELATIONSHIP BETWEEN [OIII] 
FWHM AND 

Whatever the cause of the systematic deviations in Fig. 

1, it possible to remove the systematic effect empirically 
by using a non-linear transformation from uy^Q ///jjcorr to 
cr*. For type-1 AGNs a corrected estimate of the velocity 
dispersion can be defined by the equation: 

logq[o///]]corr = (0.67 ± 0.09) log[FVKiIM[o ///] /2.35] 

-l-(0.74±0.02) (1) 

where the coefficients have been chosen so that an OLS- 
bisector fit to cr, and (J]^o iii]corr. gives the equation cr* — 
<J[0 iii]corr.- The comparison between the measured ct, 
values and the estimates using Eq. 1 are shown in Fig. 

2. The scatter for the type-1 AGNs has been reduced 
to ±0.12 dex. Type-2 AGNs cannot be used to get M, 
via the Dibai method and hence will not be considered 
further, but the corresponding transformation for type-2 
AGNs has a slope of 0.82 ±0.11. This is consistent with 
the slope of Eq. 1, but there is an offset of 0.06 dex in 
the sense that Eq. 1 would slightly over predict cr, for 
type-2 AGNs. For radio-quiet type-2 AGNs the scatter 
is reduced to ±0.12 dex, and even for radio- loud AGNs 
it is reduced to ±0.14 dex. 



o 




1.6 1.8 2.0 2.2 2.4 2.6 2.8 3.0 

Log (^-uncorr 

Fig. 3. — Th e M, -(T* relationship for the radio-quiet AGNs 
of Shi elds et al.l ([2003) (shown as crosses) and for the non-BLR 
mass estimates compiled by IGiiltekin et al.l ||2009| ) (shown as 
open circles). For the Shields et al. (200j|) AGNs c* has been 
taken simply to be FWHM[o /2.35 with no correction. The 
solid circles joined with short vertical lines t o the open cir- 
cles below them represent th e M» revisi ons of [ Humphrev et alj 
l|2008), Gebhardt & Thoma^ ||2009|) . and Ivan de n Bosch (200^ 
A typical error bar for the AGNs is shown at the upper left. 
The error in (y^uncorr is the scatter found by Nelson &; Whit tle 
lUpBI), and the error in M » has been taken to be ±0.25 fol- 
lowing Bochkarcv & Gaskelll ||2009IV The diagonal line is the 
[Gultekin et al.. 1, 2009, ) fit. 

4. THE M.-CT* RELATIONSHIP FOR AGNS 

[Shields et al.l (|2003l ) made estimates of M, and cr* for 
a large sar nple of AGNs with o bs ervations by[B rothertoj 
(ll996allbri. liSrupe et all (|1999D . iMcIntosh et al.i (|1999fK 
iDietrich et al. I (|2002[ ). and themselves. This sample 
has an adva n tage over the nearby AGN sample of 
iGreene fc Hoi (|2006[) in having AGNs covering the up- 
per end of the relationship defined by normal g alaxies 
(see IGiiltekin et al.l[2009f l. IShields et aP (|2003| ) calcu- 
lated black hole masses from FWHMs of H/3 by the Dibai 
method, and cr* was just taken to be FWHM[o/7/]/2.35. 
Details of th e treat ment of the data can be found in 
IShields et aP (|2003f ). Fig. 3 shows the M.-ct, rela- 
tionship for these BLR blackhole mass estimates and 
also for the non-BL R blackhole mass estimates com- 
piled bv IGiiltekin et al.l (j2009l) . For three galaxies re- 
vised masses are also shown (see discussion below). M, 
estimates from BLR kinematics have always been nor- 
mali zed to the st e llar-dy namical M, -cr* relationship. 
The IShields et all (l200l masses have been renormal- 
ized to the updated IGiiltekin et al.l (|2009f ) relationship. 
The renormalization has been made over the range 
2.0 < \og (J [[o II i]]corr < 2.4 whcrc the corrections to 
O HI] estima tes of cr* are smallest (see Fig. 1). 
Shields et al.l (2003) masses need to be decreased 



the 
The 

by 0.11 dex to minimize their median residual from the 
IGiiltekin et al.l (2009.) line in this range. This small 
renormalization has a negligible effect on t he results. 

It c an be seen from Fig. 3 that the IShields et al.l 
(|2003[ ) objects dev iate systematically fro m the relation- 
ship given by the IGiiltekin efall ([2009') galaxies. For 
log (J* uncorr < 2.2 the inferred masses of the.Shields et al] 
(|2003f) galaxies all lie above the lGultekin et all (|2009f ) fit. 
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Fig. 4. — The M, -a* relationship for the BLR and non-BLR 
mass estimate samples as in Fig. 3 but with correction of the stellar 
velocity dispersion estimates for the AGNs as described in the text. 
For clarity the three older mass estimates in Fig. 1 that have new 
values have been omitted. The dashed red curve is Eq. 2. The 
dashed horizontal error bar in the upper right indicates the effect 
of the uncertainty in the slope in Eq. 1. AH other symbols have 
the same meaning as in Fig. 3. The error in log cr.corr has been 
taken to be the scatter in Fig. 2. 

and for log cr^unco rr > 2.7 they a ll lie below. 

Fig. 4 plots the [Shields et all (|2003D galaxies with the 
corrected estimates of the stellar velocity dispersions 
given by Eq. 1. It can be seen that th e scatter is less 
and the systematic differences from the iGiiltekin et al.l 
HOOQ) points are reduced. The agreement is particularly 
good at the high-mass end. The typical error bar (up- 
per left corner) is now smaller. The uncertainty in the 
slope of Eq. 1 introduces an additional uncertainty in cr, 
for the highest and lowest mass black holes. The am- 
plitude of this uncertainty for a typical high-mass black 
hole is indicated by the horizontal error bar at the upper 
right. One of the free parameters in the Dibai method 
is the slope of the relationship between the effective ra- 
dius, R, of the BLR and the optical luminosity, Lagn- 
This introduces an additional uncertainty i n the - 
relationship for high- and low-mass AGNs. [Shields et all 
p003. ) took R oc LP-^ . This is in goo d agreement with th e 
best current estimate of i? cx L^-^"^ (jBentz et al.ll2009bD . 

5. DISCUSSION 

5.1. The masses of the most massive black holes 

The revised AGN M, - cr* relationship of Fig. 4 sup- 
ports recent upwards revisions of stellar dynamical A/, 
estimates in massive galaxies. Two effects are respon- 
sible for the revisions. The first is that including dark 
halos in stellar d ynamical modelling of the most mas- 
sive galaxies (see iGebhardt fc Thomas! l2009t ) increases 
the estimated masses of their black holes by about a 
factor of two. The masses of two of the most massive 
black holes have been increased already (see Fig. 3), and 
more masses of the most massive black holes are expected 
to increase similarly (K. Gebhardt - private communica- 
tion) . The second effect is that triaxial orbit calculations 
can give a higher A/, than axisymmetri c models. This 
has be en de monstrated for NGC 3379 b y Ivan den BoschI 
(|2008D and Ivan den Bosch et al.l (|2009[ ). This is likely 



to be a general effect and triaxial galaxies are believed 
to dominate among high-mass galaxies. The two sys- 
tematic effects are independent and are expected to act 
together for some high-mass galaxies. Further upward re- 
visions of the masses of the most massive black holes are 
thus likely. Some support for this already comes from di- 
rect measurement of cr* in AGNs with black hole masses 
determined from reverberation mapping. For some of 
the AGNs with loga* ~ 2.3 their masses can be up to 
~ 0.7 dex above the stellar dynamical M» relation- 
ship pasvra et al.ll2007l : I Watson et al] 120081: IWoo et all 
l2008l) . 

5.2. Curvature in the Af, -cr, relationship 

It can be s een fr om Fig. 4 that f or both the 
IGiiltekin et al.l (|2009[ ) galaxies and the IShields et al.l 
(12003D AGN^the Giiltekin et al. fit underpredicts log A/, 
by about 0.5 dex for galaxies with logAf, > 9. Even be- 
fore the recent upwards revision of mass estimat es of the 
highest mass black holes (see previous section) , IWvithd 
(i2006i) had argued that a curved log-quadratic relation- 
ship was a better fit to the Af, relationship than a 
simple power-law. iHul (|2008l ) also suggests upwards cur- 
vature of the relationship at the high-mass end. 

It is now recognized that the shape of the Af, ~ cr, re- 
lationship defends on the type of galaxy being looked at 
(jGrahaml 12001 lHull2008h . Defining the form of the re- 
lationship for a mixed sample of galaxies is particularly 
problematic at the lower end. It will take further work to 
establish the true shape of the top of the M, - relation- 
ship fo r AGNs, but if we assume that the IGiiltekin et aTl 
(j2009| ) line is a good representation for log a* < 2.5, then 
the median slope of a power law from log cr* = 2.5 to the 
top ten points in Fig. 4 is 5.5. For computational con- 
venience in interpreting observations, these two lines can 
be approximated by the log-quadratic: 

logAf. = 1. l(log cr*)2 -0.21ogcr* -f 2.82 (2) 

It should be noted that, for the points in Fig. 4, the 
reduction in per degree of freedom in going from a 
powe r-law to a log-quadrat ic relationship is not signifi- 
cant. IGiiltekin et al.l (|2009D get a similar result for their 
galaxies alone. This is because of the scatter in the lower 
half of the relationship. Whether or not a log-quadratic 
curve is formally needed depends on the choice of low 
fj* galaxies included. The significant thing about Fig. 4, 
however, is that the slope of the A^f, -ct* relationship at 
the hi gh-mass end is steeper than in th e IGiiltekin et aL 
('2009') fit and similar previous fits fe.g.. iTremaine et al. 
2002). 

Since the determination of the slope of the upper 
part of the M, - cr* relationship (see above) depends 
in part on stellar dynamical mass estimates, some of 
which are expected to increase (see section 5.2), the 
slope of the top of the relationship could well increase 
further. Even though curvature in the M,-a^ rela- 
tionship is currently only marginally statistica lly sig- 
nificant (see discussion in IGiiltekin et al.l I2009D , there 
are good reasons to believe that the relationship must 
curve upwar ds strongly. This is be cause it has long 
been known (jOegerle &: Hoesseilll991[ ) that the relation- 
ship between galaxy luminosity and cr* (|Faber fc JacksonI 
Il976f ) curves strongly at the high-luminosity end and a* 
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Fig. 5. — The dispersion logcr* (in standard deviations) as a 
function of log M, . The bins have been chosen to have equal num- 
bers of objects. 

sat urates. The highest veloc ity dispersion galaxies found 
by (jOegerle fc Hoesselll991[ ) have logcr* ~ 2.6, and fur- 
ther searches have failed to turn up any galaxies with 
log cr* > 2.65 (ISalviander et al.ll2008f ). The velocity dis- 
persions for the most massive black holes in Fig. 4 are 
consistent with this limit. The steepening in Fig. 4 re- 
moves the discrepan cy between the M. -g^. and M,- 
Lbuioe relation ships (iBernardi et al] l2007bt iLauer et al.l 
I2007f) . and as IWvithd (| 2006^ points out, the steepen- 
ing increases the local density of the highest mass black 
holes. 

5.3. Dispersion in the -cr* relationship as a 
function of M, 

Another interesting thing in Fig. 4 is the tightness 
of the A/, - (T* relationship at the high-mass end. This 
impression is confirmed in Fig. 5 which shows the dis- 
persion in loga^ for the combined samples in Fig. 4 as 
a function of logAf,. It can be seen that the disper- 
sion for the highest-mass objects is significantly lower 
than for the remaining objects. An F-test gives a prob- 
ability, p = 0.01 that the dispersions of objects with 
M, > 10^ and lO'^-^ < M, < lO^Mg are drawn from 
the same parent population. The decrease in dispersion 
at high masses is se en for the BLR estimates and the 
iGiiltekin et all l2009l estimates separately, but at lower 
significance. 

The scatter at the lower en d of the M« - g^, rel ation- 
ship is similar to that found by I Greene "fcHo (2006) from 
direct measurements of the stellar velocity dispersion by 
for AGNs with BLR masses determined by the Dibai 
method. It is also similar to that found bv lOnken et al.l 
((2004) for reverberation-mapped AGNs with a mixture 
of direct measurements of stellar velocity dispersions 
and estimate s from FWHMfo //7]/2.35). In Fig. 1 of 
I Greene fc Hd (|2006l ) the increase in scatter going from 
the high-mass galaxies with stellar dynamical determi- 
nations to the lower-mass AGNs is particularly striking. 
The results presented here argue that this effect is real 
and not a consequence of larger measuring errors for the 
AGNs. 



The increase in scatter in the M, - cr* relationship as 
one goes to lower black hole and galaxy masses proba- 
bly tells us something about blackhole growth associated 
with the different processes that dominate the evolution 
of galaxies of different masses. Elliptical galaxies and 
classical bulges in spirals are the result of violent merg- 
ers. For lower mass galaxies, secular evolution domi- 
nates instead. Secular evolution can be recognized by 
the presence of bars, pseudobulge s, and other features 
(see lKormendv fc Kennicuttll2004f ). It is now know that 
galaxies with bars and pseudo-bulges have systematically 
lower blackhole masses than predicted by the M, - a^, re- 
lationship for ellipticals and classical bulges (jHul 120081 : 
iGrahamI 120081 ). Secular evolution must therefore be a 
cause of some of the increase in scatter in Fig. 5. One 
test of whether it is the only cause of the trend is to 
look at the dispersion in logcr* for the pure ellipticals. 
Inspection of the ellipticals in the Giiltekin et al. sam- 
ple shows that the dispersion rises from » 0.05 dex at 
logAf, « 9.5 to « 0.11 at logM, « 7.5. Unfortunately, 
the significance of this trend is low (p — 0.25) because of 
the small sample size. If the dispersion does indeed rise 
for the lowest mass ellipticals, this would suggest that 
secular evolution at low masses is not the only cause of 
the dependence of dispersion on black hole mass. 

For the highest mass galaxies the low dispersion in 
M, -<T* is probably associated with the dichotomy be- 
tween elliptic als with "cores" and those with central 
"extra hght" (|Kormendv et al.ll2009D . "Core" elhpticals 
have had their cores scoured out by supermassive black 
hole binaries (Ebisuzaki et al. 1991) after "dry" dissipa- 
tionlcss mergers. "Extra-light" ellipticals have had major 
star formation in their centers following "wet" starburst 
mergers. Ellipticals w ith M, > 10^ are core ellipticals. 
IBernardi et al.l (j2007al) have shown that the brightest el- 
lipt icals have a very low scatter i n the fundamental plane, 
and' Kormendv fc Bendeil (|2009D have shown that core el- 
lipticals have a tight correlation between cr* and the stel- 
lar light deficit. All these things and the lower dispersion 
in the Af, - tr, relationship point to core ellipticals having 
very uniform properties. 

5.4. Evolution oftheM^-a^, relationship 

The radio-quiet AGNs in the lShields et al.l (120031) sam - 
ple go up to a redshift of 3.286. iShields et all (j2003f ) 
claimed that deviations from the M, - relationship 
as a function of redshift showed that the relationship 
was unchanged out to z « 3 and hence consistent with 
the growth of massive bulges and black holes occur- 
ring simultaneously. The results presented here support 
this conclusio n even though the analyzes are different. 
IShields et al.l (|2003[ ) used uncorrected estimates of cr* 
(see Fig. 3) which, as discussed above, are systematically 
too high for t he most massive AGN s. To look for evolu- 
tionary effects IShields et al.l (|2003D looked for a redshift 
dependence in residuals from the iTremaine et al.l (|2002f ) 
Af, -tJ* relationship. The slope of this r elationship is 
somewhat flatter than the IGiiltekin et al.l ((2009,) slope 
(4.02 versus 4.24) and significantly flatter than the slope 
of 5.5 for th e upper end of the re lationship in Fig. 4. The 
flatter slope IShields et all (|2003D adopted for the M,-a^ 
relationship compensated for the larger estimates of 
at high masses. 

Although the overall Af, - tr* relationship seems sim- 
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ilar at low and high redshifts much work needs to be 
done. The relationship really needs to be compared, not 
for heterogeneous collections of galaxies, but for differ- 
ent evolutionary stages of galaxies that will end up simi- 
lar. Although there does not seem to be strong evolution 
of the very highest mass galaxies and their black holes, 
strong evoluti on of lower mass g alaxies and black holes is 
quite possible ()Woo et al.ll2008D . since there is more scat- 
ter in the M, - at, relationship. In studying evolution, 
however, care needs t o be taken to allow fo r selection 
biases (see for example ISalviander et al.ir2007l ). 

6. CONCLUSIONS 

An improved transformation of FWHM[o///] to cr* 
has been proposed. This makes the velocity disper- 
sions of the highest mass AGNs consistent with the up- 
per limit on cr* in the nearby universe and requires a 
steepening of the M, - tr* relationship for the most mas- 
sive galaxies. The dispersion about the M, -cr* rela- 



tionship decreases as the mass of black holes increases. 
The tightness of the relationship for the most massive 
AGNs provides strong support both for the reliability 
of the proposed FWHM[o ///] to cr* transformation and 
of the Dibai method of estimating black hole masses 
from single-epoch spectra of type-1 AGNs. There is no 
evidence yet evolution of the M, - a* relationship from 
z ~ 3 to the present day but much further work on this 
is needed with better-defined samples. 
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